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ABSTRACT 


The  author  uses  two  equations  to  relate  velocity,  density,  flow  rate, 
and  travel  time  in  a  traffic  stream.  One  is  the  familiar  equation  of 
continuity,  the  second  is  an  integral  equation  expressing  the  conser¬ 
vation  of  fluid  in  terms  of  travel  times  in  the  fluid.  A  situation  is 
investigated  where  time -dependent  flow  rates  into  a  bottleneck  tempo¬ 
rarily  exceed  its  capacity.  Expressions  are  found  for  queue  sizes, 
the  location  and  velocity  of  shock  waves,  and  delays  to  travellers  in 
the  stream. 
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TRAVEL  TIMES  THROUGH  SHOCK  WAVES 


1.  The  Conservation  of  Fluid 

In  recent  years  Lighthill  and  Whitham  (1955),  Richards  (1956),  De  (1956), 
Bick  and  Newell  (i960)  and  Pipes  (1961)  have  discussed  the  propagation  of  shock 
waves  in  a  traffic  stream.  It  Is  assumed  that  the  vehicles  can  he  replaced  hy 
a  continuous,  compressible  fluid.  It  is  then  shown  that  the  flow  rate,  q(x,t) 
and  the  density,  k(x, t)  ,  satisfy  the  familiar  equation  of  continuity. 


(1) 
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The  dimensions  of  q(x,t)  are  vehicles  per  unit  of  time  and  the  dimensions  of 
k(x,t)  are  vehicles  per  unit  distance.  In  high-density  regions  it  is  usually 
assumed  that  the  velocity,  v(x,t)  is  only  a  function  of  the  density;  hence 
the  flow  rate, 

(2)  q(x, t)  =  k(x,t)v(x,t) 

need  only  he  written  as  a  function  of  the  density.  The  derivation  of  (l)  follows 
from  the  fundamental  law  that  the  traffic  fluid  is  conserved  at  all  points  along 
the  road  while  special  cases  of  (2)  have  been  conjectured  on  theoretical  as  well 
as  experimental  grounds.  In  the  sections  which  follow  we  assume  that  q(x,t) 
in  (2)  has  a  well-defined  hut  not  necessarily  unique  inverse. 

Solutions  of  (1)  and  ( 2 )  satisfying  appropriate  boundary  conditions  at 
time  zero  have  been  found  for  certain  flow  patterns  in  the  traffic  stream. 
Lighthill  and  Whitham  (1956)  discuss  the  qualitative  characteristics  of  the 
propagation  of  shock  waves  through  the  fluid.  Richards  (1956)  assumes  the  linear 
Greenshields  (1945)  relation  between  velocity  and  density  in  (2)  and  from  (l) 


finds  the  rate  at  which  discontinuities  in  the  stream  density  (i.e.,  shock  waves) 
propagate  away  from  an  intersection.  De  (1956)  finds  explicit  solutions  for  the 
propagation  of  shock  waves  through  a  bottleneck  whose  capacity  falls  off  linearly 
with  time,  reaches  a  minimum  value  and  then  increases  linearly  with  time  to  its 
norroa_L  value.  In  each  of  these  papers  the  authors  assume  that  the  density-flow 
reuaoions  of  (2)  which  have  been  found  to  hold  under  steady-state  flows  are  also 
valid  for  the  time -dependent  situations.  Pick  and  Uewell  (i960)  extend  these 
rGsmts  'ey  studying  tnc  j.n*ceraction  between  two  opposing  streams  of 
tiaffic  wnen  the  velocity  in  one  stream  depends  on  the  density  in  the  opposing 
as  well  as  its  own  stream. 


Fig.  I  -  THE  MOTION  OF  A  SHOCK  WAVE. 

A  lucid  derivation  of  the  formation  and  motion  of  a  shock  wave  has  been  given 
by  Pipes  (l9bl);  in  this  paper  he  also  points  out  that  the  Li gh thill -Whit ham  and 
Richards  theory  of  shock  waves  are  equivalent.  Consider  the  low  (l)  and  the 
high  (2)  density  regions  in  Figure  1.  The  interface  separating  these  two  regions 
is  a  shock  front  which  travels  with  a  velocity,  u  .  The  speed  of  the  stream 
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in  region  (l)  relative  to  the  shock  wave  is  v  -  u  ;  the  relative  speed  in 
region  (2)  is  vg  -  u  .  Since  the  net  flow  rate  of  vehicles  into  the  wave  must 
equal  the  net  flow  rate  out  of  the  wave,  we  have 

(3)  \Cvj_  -  u)  =  k2(v2  -  u)  . 

The  velocity  of  the  shock  wave  is  therefore  equal  to  the  ratio  of  the  differences 
in  flow  to  the  differences  in  density  in  each  region.  © 


v 


u  = 


2^2 


k2  ‘  kl 


Vlkl  *2  ‘  *1 


k2  "kl 


In  the  limiting  case  of  differential  arguments  where  k  and  q  are  continuously 
varying  quantities  we  have  u  =  Aq/Ak  — ►dq/dk  .  If  we  consider  the  special  case 
where  is  zero  hut  where  k£  is  large  and  k  is  negligible,  we  find  that 

the  shock  wave  moves  backwards  (negative  velocity)  with  a  velocity  equal  to  the 
incoming  flow  rate  divided  by  the  density  in  the  stalled  region.  In  other  words, 
the  inventory  of  vehicles  which  is  increasing  at  a  rate  q1  is  stored  in  a  space 
equal  to  the  number  of  vehicles  t:mes  the  average  spece,  k”1  ,  occupied  by  each 
vehicle.  The  velocity  of  the  shock  wave  simply  corresponds  to  the  rate  at  which 
inventory  space  must  be  provided  to  accomodate  vehicles.  If  the  flow  rate  can 
be  written  solely  in  terms  of  the  density.  Equation  (l)  can  be  rewritten  in  the 
form. 


(5)  "WM-o  . 

Pipes  has  shown  that  (l)  can  be  written  in  terms  of  flow  rates  by  substituting 
q  for  k  in  (5). 
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Just  as  there  is  an  equation  which  relates  densities  and  velocities  to 
flow  rates  in  a  traffic  stream  there  is  also  one  which  relates  travel  time  to 
densities  and  flow  rates,  loth  depend  upon  the  principle  of  conservation  of 
matter;  while  one  can  be  derived  from  the  other  they  do  not  always  provide 
identical  solutions. 


*  / 

Origin  Destination 


Fig.2- TRAVEL  FROM  x  TO  y 

-nc  time  to  travel  a  fixed  distance  in  the  traffic  stream  can  be  obtained 
as  follows;  Consider  Figure  2  where  two  lines,  x  and  y  ,  are  drawn  across  a 
road.  At  tjjne  t  ,  an  observer  at  y  makes  a  count  of  the  total  number  of 
vehicles  in  the  interval  (x,y)  .  The  traffic  stream  flows  from  x  to  y 
^ i . e . ,  from  xeft  to  right);  at  a  time  t  +  T(x,y,t)  the  vehicle  located  at 
x  at  time  t  passes  the  observer  at  y  ,  T(x,y,t)  is  the  travel  time  of  a 
vehicle  going  from  x  to  y  with  its  trip  originating  at  the  time  t  .  Since 


the  total  count  of  vehicles  along  the  road  at  time  t  must  equal  the  total 


flow  past  the  observer  in  the  interval  (t,  t  +  T)  ,  we  have 


(6) 


t+T(x,y,t) 


k(r,t)dr  = 


7 


q(y,s)ds 


If  the  flow  rate  is  a  constant  and  if  the  density  is  only  a  function  of  x  , 
i.e.,  stream  velocity  is  not  primarily  constrained  by  stream  density,  the 
explicit  solution  of  (6)  is  independent  of  the  arrival  time  t  and  only  depends 

upon  the  length  of  the  trip. 

* 

The  relation  between  the  integral  and  the  differential  formulation  of  the 
conservation  of  fluid  is  easy  to  establish  by  taking  the  partial  derivative  of 
Equation  (6)  with  respect  to  y  and  observing  that  the  velocity  of  the  stream 
at  a  point  y  and  a  time  t  +  T(x,y,t)  equals  the  reciprocal  of  the  partial  of 
the  travel  time  with  respect  to  the  trip  destination,  i.e., 


(7)  v(y,  t  +  T)  = 

Since  solutions  of  (5)  can  be  expressed  in  terms  of  u(k)  as  well  as  position 
and  time,  it  should  also  be  clear  that  Equation  (6)  can  be  rewritten  in  terms 
of  the  speed  of  shock  waves.  • 

2.  The  Motion  of  a  Shock  Wave 

In  this  section  we  use  Equation  (4)  to  predict  the  motion  of  a  shock  wave 
which  develops  upstream  from  a  bottleneck  and  influences  the  flow  rate  of  the 
traffic  stream  before  it  arrives  at  the  bottleneck.  In  the  following  section 
these  solutions  are  used  in  Equation  (6)  to  predict  the  travel  time  through  the 
shock  wave. 
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Fig. 3-  THE  POSITION  OFA  SHOCK  WAVE. 

Consider  the  flow  rate  into  the  bottleneck  of  Figure  3,  The  beginning  of 
the  bottleneck  is  located  at  x  =  y  ;  for  our  purposes  x  =  0  marks  the  origin 
of  trips  through  the  stream  and  also  provides  a  convenient  reference  point  where 
flow  rates  and  densities  can  be  measured.  The  stream  moves  from  left  to  right; 
the  flow  rate  at  the  origin  is  q(Q,t)  -  q(t)  and  the  density  is  k(x,t) 

=  k(t)  .  At  time  zero  we  assume  that  the  flow  rate  is  low  enough  so  that  the 
constant  capacity  flow  rate  of  the  bottleneck  is  greater  than  the  flow  rate, 
q(y.,t)  .  If  ve  neglect  statistical  variations  in  the  arrivals  of  vehicles  at 
the  bottleneck  and  if  we  assume  that  the  flow  rate  into  the  bottleneck  is  less 
than  its  capacity,  no  queue  of  vehicles  forms.  As  the  flow  rate  increases  with 
time  and  exceeds  the  capacity  of  the  bottleneck,  a  queue  forms  and  grows  back 
into  the  stream.  As  the  time -dependent  flow  rate  peaks  and  begins  to  fall  off, 
the  queue  will  continue  to  increase  until  the  flow  rate  into  the  queue  is  less ' 
than  the  capacity  flea;  rate  out  of  the  bottleneck.  The  queue  will  vanish  some 
time  after  peak  flows  have  been  observed  at  the  origin. 


Cumulative  flows  are  denoted  by  capital  letters,  i.e.. 


If  the  fluid  travelled  infinitely  fast  from  the  origin  to  the  bottleneck, 
we  would  find  that  the  size  of  the  queue,  i(t)  ,  would  equal  cumulative  net 
flows  into  the  bottleneck. 

(8)  Kt)  =  Q(t)  -  Q(sq)  -  C(t  -  sQ)  sQ  £  t 

where  sQ  is  the  instant  when  flow  rates  first,  exceed  C  .  If  we  also  assume 
that  each  vehicle  occupied  a  space  in  the  queue  equal  to  the  reciprocal  of  the 
jam  density,  k^  ,  the  position  of  the  shock  wave  would  represent  the  head  of 
the  queue  and  would  be  located  at 

(9)  x(t)  =  y  -  k"1  i(t)  sQ  1  t 

The  smallest  value  of  x(t)  would  occur  at  the  time  when  the  queue  size  has 
its  largest  value,  i.e.,  that  instant  of  time  t  =  s1  which  is  the  smallest 
root  of 

(10)  q(t)  -  C  =  0 

On  the  other  hand,  if  we  include  the  effects  of  finite  stream  velocity, 
the  flow  rate  at  the  bottleneck  first  equals  the  capacity,  C  ,  at  a  time 
tQ  >  sQ  when 

(n)  q.(y,tQ)  =  C 

If  we  further  assume  that  outside  the  queue  the  stream  velocity  is  the  free 
velocity,  vf  ,  and  is  independent  of  position,  the  flow  rate  at  a  point  x  is 

C12)  l(x,t)  =  q(0,t  -  v^x)  =  q(t  -  v^x) 
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In  other  words,  the  flow  rate  function  at  x  is  equal  to  the  flow  rate  function 
at  the  origin  shifted  "by  the  time  to  travel  from  0  to  x  .  The  solution  of 
(ll)  becomes, 

U3)  =  s0  +  ^f1  * 

Consider  that  there  are  only  two  density  regions:  one  where  vehicles  are  outside 
the  queue,  the  other  where  vehicles  are  closely  packed  in  the  queue  that  forms 
upstream  from  the  bottleneck: 


m 

(14a) 

k(x,t)  =  v”Xq(x,t ) 

x  outside  the  queue 

(14b) 

=  kJ 

x  inside  the  queue 

Substituting  (l4)  for  k^  and  k^  and  q(x,t)  and  C  for  q1  and  q^ 
in  (4),  we  obtain  the  velocity  of  the  shock  wave,  in  terms  of  its  co-ordinate  at 
time  t  . 


(15) 


/ .  \  dx 
U(t)  "  dt 


q.(t  -  v^x)  -  c 

v^Xq(t  -  v"Xx)  -  kj 


This  first  order  differential  equation  is  nonlinear  in  x  and  its  derivative. 
While  numerical  solutions  can  be  found,  a  useful  analytic  expression  for  the 
position  of  the  shock  wave  as  a  function  of  time  can  be  obtained  by  expanding 
the  right-hand  side  of  (15 )  about  t  and  neglecting  terms  smaller  than  xv~2  . 
The  result  is  a  first  order  linear  differential  equation  of  the  form. 


(16) 


-ff  +  a(t)x(t)  =  b(t) 
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where  a(t)  ,  b(t)  are  time -varying  coefficients  involving  the  flow  rate  and 
its  derivatives.  The  complete  solution  of  (l6)  satisfying  the  initial  condition 
that  the  shock  wave  he  located  at  x  =  y  when  the  flow  rate  at  the  bottleneck 
first  exceeds  C  is 


(17) 


x(t)  = 


C(t  -  t_)  -  (Q(t) 


k.  -  k(t) 

J 


Q(t0))  k  -  k(t  ) 

- -  +  -J -  y 

kj  -  k(t) 


where  k(t)  -  k(0,t)  and  t^  is  given  by  (13).  The  shock  wave  reaches  its 
extreme  upstream  point  at  time  t^  when 


(l8)  x(tx)  =  ( <1  ( t;L )  -  C)vf  _1 

and  returns  to  the  bottleneck  at  t^  when  x(t^ )  —  y  .  By  substituting  y  for 
the  left-hand  side  of  (l4)  we  find  that  t  is  explicitly  independent  of  the 
jam  density, 

te)  Q(V  -  Q(t0)  "  c^2  "  to)  =  -  t(t0)] 

By  neglecting  terms  smaller  than  yvf2  we  find  that  this  return  time,  t  ,  can 
also  be  obtained  by  equating  the  cumulative  flows  into  the  bottleneck  with  the 
maximum  flow  rates  which  can  pass  through  the  bottleneck  in  the  capacity  con¬ 
strained  interval  (tQ  ,  t  )  ,  i.e., 

(20)  Q(t2  -  yv"1)  -  Q(t0  -  yv”1)  =  C(t£  -  tQ)  . 

In  other  words  the  duration  of  the  capacity  constrained  interval  is  not  (as  a 
first  approximation)  affected  by  the  finite  size  of  the  vehicles  or  by  the  average 
speed  of  the  vehicles  within  the  queue  itself  even  through  the  shock  wave 
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influences  upstream  densities  and  flew  rates  during  the  course  of  its  progress. 

If  the  density  of  vehicles  at  the  origin 'at  time  t  is  small  and  close  to  the 

density  at  time  tQ  ,  the  second  term  in  (14)  is  approximately  equal  to  y  ; 

the  first  term  is  approximately  equal  to  kj1^)  -  Q(tQ)  -  C(t  -  tQ)]  which, 

except  for  the  factor  k^1  ,  is  identical  to  Equation  (8). 

Figure  4  is  a  plot  of  the  position  of  the  shock  wave  which  forms  upstream 

from  the  "bottleneck  in  Figure  3.  The  flow  rate  at  the  origin  is  q(t) 

2  -*t 

=  3,000  t  e  cars/hour  where  t  is  measured  in  hours.  The  bottleneck  is 
five  miles  downstream  from  the  origin,  C  -  1450  cars/hour  and  k  =  220  cars/mile. 
The  solid  lines  are  plots  of  the  shock  wave  for  several  values  of  the  free 
velocity  vf  and  its  distance  from  the  origin  in  miles  is  shown  on  the  right- 
hand  ordinate;  the  dashed  line  is  a  plot  of  the  flow  rate  q(t)  with  values 
given  on  the  left-hand  ordinate. 

As  vf  *  00  the  position  of  the  shock  wave  is  given  by  Equation  (9);  it 
starts  at  tQ  =  0.700  hour  and  has  its  smallest  value  when  t  =  1.00  hour.  As 
the  velocity  of  the  stream  outside  the  queue  is  changed,  three  effects  can  be 
observed:  (i)  The  shock  wave  starts  later  as  vf  is  reduced  because  it  takes  a 
longer  time  for  vehicles  to  travel  from  the  origin  to  the  bottleneck,  (ii)  the 
largest  upstream  position  of  the  shock  wave  increases  as  v  gets  smaller 
(iff)  the  width  of  the  peak  narrows  as  v^,  decreases. 

3.  Travel  Time  Through  the  Shock  Wave 

Consider  the  time  to  make  a  trip  from  the  origin  to  the  bottleneck  in 
£xgure  3*  Part  of  the  trip  is  made  in  the  low-density  region,  part  in  the 
high-density  region,  i.e., 

(total  travel  time)  =  (time  to  travel  to  shock  wave) 

+  (time  to  travel  beyond  shock  wave) 
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POSITION  OF  THE  SHOCK  WAVE 


Substituting  (l4)  into  (6)  for  t  >  t  Ve  get 


(21)  =  c  1 J  v^qft  -  uv^du  +  f  k  du  "j- 

0  xft)  J 

=  C"1  |  Q(t)  -  Q(t  -  v"1x(t ) )  +  kj (y  -  x( t) )  j-  . 

Equation  (21)  can  be  written  in  two  alternative  forms.  Expanding  about 
t  ,  neglecting  terms  0(vf  )  ve  can  express  the  travel  time  in  terms  of  the 
normalized  density  rj(t)  =  k(t)/k.  £  1  , 

J 

T(°>y»t)  =  C  H (y  -  x(t)  +  T)(t)x(t)J  t  ^  t 

In  other  words,  the  total  travel  time  is  the  time  to  travel  through  a  queue 
moving  at  the  constant  velocity  plus  a  term  which  is  proportional  to 

the  product  of  the  nromalized  density  and  the  position  of  the  shock  wave  at  the 
origin  of  the  trip. 

From  (21)  we  also  have 

(23)  T(0,y,t)  =  C-1y  a(t)k(t)  +  (l  -  a(t))k. 

J 

where  a(t)  is  the  fraction  of  the  trip  covered  by  the  high-density  or  slow 
moving  region  at  time  t  .  Since  a(t)  ^  1  ,  the  term  in  brackets  is  a  density 
lying  between  the  density  observed  at  the  beginning  of  the  trip  and  the  jam 
density.  Hence  upper  and  lower  bounds  can  be  placed  on  the  total  travel  time. 

As  the  velocity  in  the  low-density  region  becomes  very  large,  i.e.,  as  v  _ ►  <x 

we  have 

(24)  T(0,y,t)  —  . 
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This  simply  states  that  the  total  travel  time  approaches  the  delay  of  a  vehicle 
joining  a  queue  of  size  l(t)  which  is  then  processed  at  a  constant  service 
rate  C  . 

We  have  used  vf  to  denote  the  ("free")  speed  of  the  traffic  stream  in 
.  the  low-density  region.  The  results  can  he  generalized  to  include  cases  where 
the  velocity  in  the  low-density  region  is  constant  hut  less  than  the  free  speed. 
For  example,  the  stream  speed  may  he  given  hy  expressions  of  the  form 
v  =  v(k)  <  vf  . 

In  generalizing  these  results  at  least  one  important  effect  would  have 


to  he  considered.  The  first  of  these  would  he  a  more  realistic  assumption  about 

the  behavior  of  the  stream  density  in  the  neighborhood  of  the  beginning  of  the 

queue.  Vehicles  will  not  immediately  slow  down  from  v  to  Ck”1  on  passing 

i  J 

through  the  shock  wave.  Instead,  they  will  slow  down  more  gradually.  Since 
the  density  will  also  change  more  gradually,  the  leading  edge  of  the  shock 
wave  will  give  way  to  a  region  where  the  changes  In  flow  with  respect  to 
density  is  large  hut  not  infinite.  Even  the  simplest  assumptions  about  vehicle 


deceleration  characteristics  seem  to  lead  to  much  more  complicated  expressions 
than  those  we  have  already  obtained;  f  or  this  reason  they  are  not  presented 
here. 
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